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Abstract 

The photoproduction of dijet events, where the two jets with the highest trans- 
verse energy are separated by a large gap in pseudorapidity, have been studied 
with the ZEUS detector using an integrated luminosity of 39 pb" 1 . Rapidity- 
gap events are defined in terms of the energy flow between the jets, such that 
the total summed transverse energy in this region is less than some value E^ VT . 
The data show a clear excess over the predictions of standard photoproduction 
models. This is interpreted as evidence for a strongly interacting exchange of 
a color-singlet object. Monte Carlo models which include such a color-singlet 
exchange are able to describe the data. 
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1 Introduction 



The production of events in hadronic collisions with two high transverse energy jets in the 
final state separated by a large rapidity interval provides an ideal environment to study 
the interplay between soft (non-perturbative) and hard (perturbative) QCD. 

The dominant mechanism for the production of jets with high transverse energy in 
hadronic collisions is a hard interaction between partons in the incoming hadrons via 
a quark or gluon propagator. The exchange of color quantum numbers generally gives 
rise to jets in the final state that are color connected to each other and to the remnants 
of the incoming hadrons. This leads to energy flow populating the pseudorapiditjQ re- 
gion both between the jets and the hadronic remnants, and between the jets themselves. 
The fraction of events with little or no hadronic activity between the jets is expected to 
be exponentially suppressed as the rapidity interval between the jets increases. A non- 
exponentially suppressed fraction of such events would therefore be a signature of the 
exchange of a color-singlet (CS) object. 

The high transverse energy of the jets provides a perturbative hard scale at each end of 
the CS exchange, so that the cross section should be calculable in perturbative QCD [1]. 

Previous studies of jets with rapidity gaps have been made in pp collisions at the Tevatron 
[2,3] and in photoproduction at HERA [4,5], where a quasi- real photon from the incoming 
positron interacts with the proton. Comparison with different Monte Carlo (MC) models 
suggested that some contribution of a strong CS exchange is required to describe the data, 
although the uncertainty on the contribution from standard QCD processes was large. 

In the analysis presented in this paper, photoproduction of dijet events with a large 
rapidity gap between jets is used to investigate the dynamics of color singlet exchange. 
The results are based on a larger data sample, than in the previous publications [4,5]. 
The MC models were tuned to better describe the data sample at the detector level. The 
CS contribution is studied and compared to MC models as a function of several kinematic 
variables and to a recent QCD-resummed calculation [6-8]. 



2 Experimental set-up 

The results presented in this paper correspond to 38.6 ±1.6 pb" 1 of data taken with 
the ZEUS detector during the 1996-1997 HERA running period. Positrons of 27.5 GeV 
collided with protons of 820 GeV, giving a center-of-mass energy of y/s = 300 GeV. 

1 The pscudorapidity 77 = — ln[tan(0/2)], where 8 is a polar angle. 
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A detailed description of the ZEUS detector can be found elsewhere [9,10]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are measured in the central tracking detector (CTD) [11], which oper- 
ates in a magnetic field of 1.43 T provided by a thin super-conducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organized in nine super-layers covering the 
polar-angl^§ region 15° < 6 < 164°. The transverse momentum resolution for full-length 
tracks can be parameterized as ct(pt)/pt = 0.0058pr © 0.0065 © 0.0014/pr, with pt in 
GeV. The tracking system was used to measure the interaction vertex with a typical 
resolution along (transverse to) the beam direction of 0.4 (0.1) cm and also to cross-check 
the energy scale of the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [12] covers 99.7% of the total 
solid angle and consists of three parts: the forward (FCAL), the barrel (BCAL) and 
the rear (RCAL) calorimeters. Each part is subdivided transversely into towers and 
longitudinally into one electromagnetic section and either one (in RCAL) or two (in BCAL 
and FCAL) hadronic sections. The smallest subdivision of the calorimeter is called a 
cell. Under test-beam conditions, the CAL single-particle relative energy resolutions were 
a(E)/E = 0.18/ VE for electrons and a{E)/E = 0.35 /VE for hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e^p. The 
resulting small angle energetic photons were measured by the luminosity monitor [13], a 
lead-scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 

3 Kinematics and event selection 

A three-level trigger system was used to select events online [10,14]. In the third-level trig- 
ger, jets were required to have a transverse energy of E 3 ^ 1 > 4 GeV and a pseudorapidity 
of r] )et < 2.5 in the laboratory frame. 

The 7p center-of-mass energy, W, and the inelasticity, y = W 2 / s, were reconstructed using 
the Jacquet-Blondel (JB) [15] method. The hadronic system was reconstructed using 
Energy Flow Objects (EFOs), which were formed by combining information from energy 
clusters reconstructed in the CAL and charged tracks reconstructed in the CTD. The 
electron (e) [16] reconstruction method was also used, in order to remove deep inelastic 
scattering (DIS) events. 

The photoproduction sample was selected by applying the following offline cuts: 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. 
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• the longitudinal position of the reconstructed vertex was required to be in the range 
-40 cm < Z VTX < 40 cm; 

• events with a scattered positron in the CAL having y e < 0.85 and E' e > 5 GeV, 
where E' e is the energy of the scattered positron, were rejected. This cut reduced 
contamination from neutral current DIS events, since the efficiency for the detection 
of the scattered positron in this region approached 100%; 

• events were required to have 0.2 < i/jb < 0.75. The upper cut on ?/jb further reduced 
contamination from the neutral current DIS events that were not removed by the cut 
on y e and the lower cut removed beam-gas events; 

• in order to reduce contributions from charged current events and cosmic-ray showers, 
events were required to have a relative transverse momentum P™ 1SS / \JE^ < 2.0 GeV 1 / 2 , 
where P™ ss and are the total event missing momentum and transverse energy, 
respectively. 

The cuts on y e and ?/jb reduced the contribution of DIS events to less than 0.5%, confined 
the phase-space region of the analysis to 0.2 < y < 0.75 and restricted the photon 
virtuality to a range of Q 2 < 1 GeV 2 with a median value of Q 2 ~ 10~ 3 GeV 2 [17]. 

Jets were reconstructed from the EFOs using the kj- algorithm [18] in the longitudinally 
invariant inclusive mode [19], which implies that any particle is included in one of the jets, 
and ordered in E 3 ^ , such that jetl had the highest E^. Events in which jets satisfied the 
following criteria were then selected: 

• jet transverse energy corresponding to E^ n > 6 GeV and E^ t2 > 5 GeV at the hadron 
level, after taking in account energy loss in inactive material and other detector effects; 

• —2.4 < rf Qt1,2 < 2.4, where rf etl and rf et2 are the pseudorapidities of the corresponding 
jets, to ensure that the jets were well reconstructed in the detector; 

• 2.5 < Arj < 4, where At] = \rf etl — rf et2 \ is the absolute difference in pseudorapidity 
between the jets; 

• 1 1 (V ctl + rf ct2 } | < 0.75, where this condition, together with the previous one, con- 
strained the jets to lie within the kinematic region where the detector and event 
simulation are well understood. 

The transverse energy in the gap, E^ AP , was calculated by summing up the transverse 
energy of all jets, without any cut on E 3 ^ , lying in the pseudorapidity region between the 
two highest-E^* jets satisfying the above requirements [20]. Gap events were defined as 
those in which P^ AP was less than an E^ VT value. The E^ VT values used in this analysis 
were E^ VT = 0.5, 1.0, 1.5, and 2.0 GeV. The gap fraction, /, was defined as the ratio of 
the cross section for gap events to the cross section for inclusive events, which pass all of 
the above cuts but have no restriction on the E^ AF value. 
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In addition, the fraction of the photon momentum participating in the hard interaction 
was calculated as x° BS = (E^ n e~^ etl + _E£ t2 e~^° t2 ) /2yE e , where E e is the energy of the 
positron beam. 

4 QCD models and event simulation 
4.1 Monte Carlo models 

The Pythia 6.1 [21] and Herwig 6.1 [22] MC generators were used to correct the data 
to the hadron level and for model comparisons. Both MCs are based on the leading order 
(LO) (2 — > 2) matrix elements together with a parton-shower simulation of additional 
QCD radiation and hadronisation models. The detector simulation was performed with 
the Geant 3.13 program [23]. 

In photoproduction interactions at LO, the photon can either participate directly in the 
hard sub-process (direct photoproduction) or first fluctuate into a hadronic state which 
then interacts via a partonic constituent carrying some fraction, cc 7 , of the photon momen- 
tum (resolved photoproduction). At leading order, therefore, CS exchange between jets 
may take place only in resolved photoproduction. For this analysis the direct, resolved, 
and CS exchange MC samples were generated separately. 

The simulation of multi-parton interactions (MPI) was included in Pythia using the 
so-called "simple mode" [21] and in Herwig by interfacing to the Jimmy library [24]. 
The minimum transverse momenta, p™ m , of the outgoing partons in the hard interaction 
and partons participating in MPI are separately adjustable in Pythia, while in Herwig 
the same parameter was used to adjust both momenta. The starting parameters for 
the tuning were taken from global fits of JetWeb [25]. The p™ in was tuned [26,27] 
for both MC programs by comparing to the data sample after the kinematic cuts were 
applied (see Section [2]). The best fit resulted in p™ m values of p™ ml = 1.9 GeV and 
p™ in2 = 1.7 GeV for Pythia and pf n = 2.7 GeV for Herwig. For both MC models the 
CTEQ5L parametrisation [28] for the proton and the SaS-G 2D parametrisation [29] for 
the photon PDFs were used. Hadronisation in Herwig is simulated using the cluster 
model [30] while Pythia uses the Lund string model [31]. 

The CS exchange is implemented in Herwig using the LLA BFKL model by Mueller and 
Tang [32]. The hard-Pomeron intercept, l + u>o, is related to the strong coupling, a s , used 
in the BFKL parton evolution by lu = a \ A [4 In (2)]. In this analysis, the default value 
of ujq = 0.3 was used. 

Pythia does not contain a simulation of strongly interacting CS exchange in hard in- 
teractions. However, a similar topology can be simulated by high-t photon exchange for 
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quark-quark scattering in LO resolved processes. Such an exchange is not expected to 
represent the mechanism of strongly-interacting CS exchange and is only used to compare 
the data to an alternative CS model. 

4.2 Resummed calculation 

The gap definition in terms of the energy flow between jets, being infrared safe, allows 
pQCD calculations to be applied. These calculations involve the resummation of large 
logarithms of Ap / . There are several sources of these large logarithms. The primary 
leading logarithms arise from soft gluon emission directly into the gap, whereas secondary 
(non-global) leading logarithms are due to emission into the gap from a coherent ensemble 
of gluons outside the gap region [7,8,33-35]. 

The calculation [6] used in this paper provides a prediction of the gap fraction with 
primary emission resummed to all orders and a correction applied for the effect of the 
clustering algorithm, and the non-global logarithms correct in the limit of large number 
of colors. The theoretical uncertainty in this calculation is estimated from varying the 
renormalisation scale between E^/2 and 2Et, where E^ is the transverse energy of the 
hardest jet. 

5 Data correction and systematic uncertainties 

The data were corrected to the hadron level, bin-by-bin, using correction factors obtained 
from a combination of direct, resolved, and CS MC samples as described in detail elsewhere 
[26,27]. 

The admixture of direct and resolved MC used in the unfolding was determined by the 
best fit to the x® BS data distribution. The combination of direct and resolved MC formed 
the non-color-singlet (NCS) sample. 

The relative amounts of NCS and CS MC used in the unfolding were determined by the 
best fit to the total energy in the gap for events in which E^ AP < 1.5 GeV, after the 
normalisation of the NCS sample was fixed using data at E^ AP > 1.5 GeV. Fitting to the 
total number of jets in the gap for events in which Etf AF < 1.5 GeV and to dcr/dE^ AF 
gave similar results. 

To correct the data the average correction factor of Pythia and Herwig was used. One 
half of the difference between those two models predictions, about 5%, was assigned to 
the systematic uncertainties. 
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A detailed study of the sources contributing to the systematic uncertainties of the mea- 
surements was performed using Herwig. The analysis cuts were varied by their respective 
resolutions estimated using Monte Carlo. 

The variation of the cuts on E^ AP and Et caused the largest contributions to the system- 
atic uncertainty. Depending upon the variable measured, their contribution ranged from 
a few to approximately 30% in regions where the statistical significance was low. 

The amount of CS exchange MC used in the unfolding was varied by ±25%, resulting in 
a variation in the cross section at the one percent level. All the above systematics were 
added in quadrature in order to calculate the total systematic uncertainty. 

The calorimeter energy scale was varied by ±3%. This uncertainty was not combined 
with the other systematics, but instead shown separately as a shaded band in the figures. 



6 Results 

The inclusive dijet cross section as a function of E^ AF is presented in Fig. [T] and Table [TJ 
At low E?f AF values, where the CS contribution should be most pronounced, the data 
demonstrate a clear excess over the NCS MC predictions. In order to estimate the amount 
of CS contribution, the direct and resolved components of each MC were mixed according 
to their predicted MC cross sections to give the NCS MC sample. The NCS and CS MC 
samples were then fitted to the data according to 

da da NCS da cs 

dE$ AP ~ 1 dE% AP + 2 dE^ AP ' 

where Pi and P 2 were the free parameters of the fit. The best fit to the data resulted in 
Pi = 1.31 ±0.01 and P 2 = 327±20 for Pythia and Pi = 1.93±0.01 and P 2 = 1.02±0.13 
for Herwig. These scaling parameters were used in this analysis when comparing the 
data to the MC predictions. The large value of P 2 for Pythia reflects the very low cross 
section of the high-t photon exchange, which is not expected to represent the mechanism 
of strongly-interacting CS exchange. The color singlet contribution to the total cross 
section, estimated by integrating the MC predictions over the entire P^ AP range, was 
(2.75 ± 0.10)% for Pythia and (2.04 ± 0.25)% for Herwig, where the errors represent 
only the statistical uncertainties of the fit. 

The inclusive dijet cross section, the gap cross section, and the gap fraction as a function 
of the separation of the two leading jets, Arj, are presented in Fig. [2] for E^ VT = 1 GeV. 
Both cross sections and gap fractions decrease as a function of Arj. In the inclusive cross 
section, both MC models with and without CS exchange describe the data equally well. 
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For the gap cross section the MC models without CS exchange fall below the data, while 
the MC models with CS exchange agree with the data. The contribution of CS exchange 
to the total gap fraction increases as the dijet separation increases from 2.5 to 4 units in 
pseudorapidity. 



CUT 
T 



Figure [3] shows the gap fraction as a function of Arj for the four values of E, 
0.5,1.0,1.5 and 2 GeV. The corresponding values are listed in Table [2j The data first 
fall and then level out as Arj increases for all values of E^ VT , although for E^ VT = 0.5 
the data are consistent with a flat distribution in At/. The predictions of Pythia and 
Herwig without CS exchange lie below the data over the entire Arj range. With the 
addition of the CS contribution, both MC models describe the data well. 

The previously published ZEUS results [4] used a different definition of the rapidity gap 
and so cannot be directly compared. The present results agree with the previous HI 
measurement [5], where the gap definition used the transverse energy in the gap as for 
the current analysis, but with slightly different kinematic cuts. The comparison is shown 
in Fig. HI where the HI data have been scaled bin-by-bin with multiplicative factors 
estimated using the Herwig MC predictions for the gap fractions at the hadron level to 
account for the difference in the phase space between the ZEUS and HI analyses. 

Figure [5] shows the gap fraction for four different values of E^ VT compared to the re- 
summed calculation [6]. The shape of the data as a function of Arj is reasonably well 
described for all values of E^ VT but the predictions lie above the data, almost everywhere 
outside of the range defined by the theoretical uncertainties. 

For comparison with other experiments and pp measurements, which are expected to 
be similar to the resolved-photon process, the cross sections and gap fraction were also 
measured as function of x° BS . These results are presented in Figs. El Hand Table[3Jfor four 
different values of E^ VT . The gap fraction decreases with decreasing x® BS and the data 
are reasonably described by both MC models only after including the CS contribution, 
especially in the resolved photon region, x° BS < 0.75, and at low E^ AF . 

The W dependence, which is important for comparison with experiments at different 
energies, is presented for the cross sections and gap fractions in Figs. [SJ, M and Table [H 
The gap fraction falls with increasing W. Both the cross sections and the gap fractions 
are described by the MC with CS included. 

The Arj and W dependencies were investigated in the resolved enhanced region. Figure [TO] 
shows the cross sections as a function of Arj in the resolved photon region, x® BS < 0.75, 
for Etf AP < 1 GeV. The gap fraction as a function of Arj is reasonably well described 
by MC models after including the CS contribution. Figure [IT] and Table [5] show the 
gap fractions as a function of Ar] for the resolved enhanced sample for the four E^ VT 
values. For E^ AP < 0.5 GeV and Etf AP < 1.0 GeV, both MC models predict almost no 
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contribution to the gap fractions from the NCS component at high values of Arj. The W 
behavior in the resolved enhanced sample is presented in Figs. [12] and [13] and Table [6] 

Although the gap fraction was measured with small errors, the difference in the model 
predictions precludes a model-independent determination of the CS contribution. 

7 Summary 

Dijet photoproduction has been measured for configurations in which the two jets with 
highest transverse energy are separated by a large rapidity gap. The fraction of events 
with very little transverse energy between the jets is inconsistent with the predictions of 
standard photoproduction MC models. The same models with the inclusion of a color- 
singlet exchange sample at the level of 2 — 3% are able to describe the data, including the 
gap-fraction dependency on E^ AP , W, x® BS and Arj. 

The difference in the model predictions precludes an accurate determination of the color- 
singlet contribution and its behavior as a function of different kinematic variables such as 
x obs or w 
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E^ AP bin ( GeV) 


cr( nb/ GeV) ± stat ± sys ± cal 


0.0-0.5 
0.5-1.5 
1.5-3.5 
3.5-7.0 
7.0 - 12.0 


1 67 ± 004 +0 - 014 +0 - 002 

V.IUI ZC U.ULKt _ .oi4 -0.006 

n i c;o + n 002 + ' 006 +0 - 000 

U.lOO ZC U.UUZ _o.006 -0.001 

210 ± 002 +0 - 009 +0 - 001 

U.Z1U ZZ U.UUZ _ 000 8 -0.002 

0.177 ±0.001 ±°oZ 
0.080 ±0.001 t° ™ t° Z 



Table 1: The measured differential cross section da/dE^ AP unfolded with the aver- 
age correction factors of Pythia and Herwig for the inclusive sample of events. 
The statistical error, systematic errors, and calorimeter energy- scale uncertainty 
on the measurement are also listed. 



Arj bin 


£CUT GeV 


/ 


± stat 


± sys 


± cal 


2.5, 


2.8 




0.053 


± 0.002 


+0.007 
-0.004 


+0.003 
-0.003 


2.8. 


3.1 


0.5 


0.047 


± 0.002 


+0.006 
-0.007 


+0.004 
-0.003 


3.1. 


3.5 


0.040 


± 0.003 


+0.008 
-0.009 


+0.002 
-0.005 


3.5. 


4.0 




0.038 


± 0.005 


+0.012 
-0.012 


+0.001 
-0.000 


2.5. 


2.8 




0.101 


± 0.002 


+0.006 
-0.005 


+0.004 
-0.005 


2.8. 


3.1 


1.0 


0.080 


± 0.003 


+0.007 
-0.005 


+0.005 
-0.004 


3.1. 


3.5 


0.061 


± 0.003 


+0.006 
-0.006 


+0.001 
-0.004 


3.5. 


4.0 




0.055 


± 0.005 


+0.014 
-0.016 


+0.003 
-0.002 


2.5. 


2.8 




0.163 


± 0.003 


+0.007 
-0.009 


+0.008 
-0.007 


2.8. 


3.1 


1.5 


0.127 


± 0.003 


+0.005 
-0.005 


+0.007 
-0.007 


3.1. 


3.5 


0.094 


± 0.003 


+0.007 
-0.005 


+0.003 
-0.005 


3.5. 


4.0 




0.092 


± 0.007 


+0.019 
-0.030 


+0.003 
-0.004 


2.5. 


2.8 




0.228 


± 0.003 


+0.011 
-0.010 


+0.012 
-0.011 


2.8. 


3.1 


2.0 


0.178 


± 0.004 


+0.012 
-0.006 


+0.010 
-0.008 


3.1. 


3.5 


0.135 


± 0.004 


+0.014 
-0.010 


+0.006 
-0.006 


3.5. 


4.0 




0.138 


± 0.008 


+0.019 
-0.035 


+0.001 
-0.009 



Table 2: The measured gap fraction f (At]) unfolded with the average correc- 
tion factors of Pythia and Herwig. The statistical error, systematic errors, and 
calorimeter energy-scale uncertainty on the measurement are also listed. 
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x° BS bin 


££UT GeV 


/ ± stat ± sys ± cal 


0.00,0.50 
0.50,0.75 
0.75,0.90 
0.90,1.00 


0.5 


o.oi7 ±0.002 j*™ ±g:Sg? 

0.018 ±0.001 t°-Z 

o.o39 ±0.002 ±8:8S 

0.272 ±0.010 t°o°o% Iffi 


0.00,0.50 
0.50,0.75 
0.75,0.90 
0.90,1.00 


1.0 


o.o28 ±0.003 ±85SS 18:88? 

029 ± 001 + 004 +0 ' 001 

U.UZJ HZ U.UUI _o.003 -0.002 

n n 7 q 4- n 002 + 005 + - 003 

U.U/J =C U.UUZ _ .005 -0.005 

4S4 ± 1 2 + - 024 +0 ' 008 

U.^d^t 3= U.U1Z _ .026 -0.008 


0.00,0.50 
0.50,0.75 
0.75,0.90 
0.90,1.00 


1.5 


047 ± 3 +0 - 005 +0 ' 001 

U.U^W =C U.UUO _ .007 -0.002 

046 ± 001 +0 006 +0 ' 003 

"• J ™ 31 U.UUI _o.005 -0.003 

14^ ± 03 +0 - 007 +0 ' 006 

U.l^O 3= U.UUO _ .010 -0.008 

o.63o ±0.015 ±8:811 18:85? 


0.00,0.50 
0.50,0.75 
0.75,0.90 
0.90,1.00 


2.0 


0.069 ±0.004 :H?o 7 18Z 

070 ± 002 +a008 +0 ' 004 

U.U/U 3= U.UUZ _ .005 -0.005 

o.227 ±0.004 ±8;8il ±8:808 
o.763 ±0.018 ±°;°i? ±8:881 



Table 3: The measured gap fraction /(x° BS ) unfolded with the average correc- 
tion factors o/Pythia and Herwig. The statistical error, systematic errors, and 
calorimeter energy-scale uncertainty on the measurement are also listed. 
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W bin ( GeV) 


££UT GeV 


/ ± stat ± sys ± cal 


150.0,180.0 
180.0,210.0 
210.0,240.0 
240.0,260.0 


0.5 


o.o77 ±0.007 t° °oll ±8:SSS 

0.049 ± 0.003 t° Z ±8:gg 
0.039 ± 0.002 t° ™ ±g : ggg 

o.o38 ±0.002 ±g5gs ±g : ggg 


150.0,180.0 
180.0,210.0 
210.0,240.0 
240.0,260.0 


1.0 


o.i45 ±0.008 j£E§ ±8:ggS 
o.o96 ±0.004 %™ ±8:SSi 
o.o69 ±0.002 j™°S ±8:SS 

062 ± 02 +0 - 006 +0 - 005 

U.UOZ 3Z U.UUZ _o.004 -0.003 


150.0,180.0 
180.0,210.0 
210.0,240.0 
240.0,260.0 


1.5 


241 ± 1 +0 - 025 +0 ' 003 

U.Z^l 3= U.UIU _o.019 -0.015 

n i c;o + n 004 +aol ° +a008 

U.IOO 3= U.UU4: _o.010 -0.004 

o.iis ±0.003 ±g58i ±8:ggg 
o.o97 ±0.003 %™ ±8:8S 


150.0,180.0 
180.0,210.0 
210.0,240.0 
240.0,260.0 


2.0 


0.338 ±0.012 i°:8§? ^ 

0.218 ±0.005 is^g ±8:8S2 

0.163 ± 0.003 iffi ±g : gg? 

o.i39 ±0.003 i°:8oi ±8:881 



Table 4: TTie measured gap fraction f (W) unfolded with the average correc- 
tion factors o/Pythia and HERWIG. The statistical error, systematic errors, and 
calorimeter energy-scale uncertainty on the measurement are also listed. 
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Arj bin 




/ ± stat ± sys ± cal 


2.5,2.8 
2.8,3.1 
o.L, 6.0 
3.5,4.0 


0.5 


021 ± 2 +0 - 003 +0 - 001 

U.UZ1 31 U.UUZ _ .003 -0.001 

nm4 + n 002 +0 - 005 +0 - 001 

U.Ul^t 3= U.UUZ _o.004 -0.001 

nmc; _i_ n nno +0.004 +0.000 
0.015 ± 0.002 _ 005 _ 0003 

0.009 ±0.003 ±8:88g 


2.5,2.8 
2.8,3.1 
3.1,3.5 
3.5,4.0 


1.0 


o.o38 ±0.002 ±8:gg2 ±g;gS 

0.024 ± 0.002 ±g;gg| j™" 
0.019 ± 0.002 +g;g« ^ 
01 fi + 04 +0 - 005 +0 - 000 

U.UIO IC U.UU^ _o.008 -0.002 


2.5,2.8 
2.8,3.1 
3.1,3.5 
3.5,4.0 


1.5 


060 ± 2 +0 - 006 +0 - 003 

U.UUU 33 U.UUZ _ .006 -0.002 

040 ± 2 +0 - 005 +0 - 003 

U.U^U 33 U.UUZ _ .005 -0.003 

027 ± 2 +0 - 006 +0 - 001 

U.UZ/ 33 U.UUZ _ .003 -0.003 

0.026 ± 0.006 j™» ±£SS 


2.5,2.8 
2.8,3.1 
3.1,3.5 
3.5,4.0 


2.0 


o.o9o ±0.003 ±g:SS ±8:33? 

n Q63 ± 003 +a008 +a003 

U.UUO 3= U.UUO _ .007 -0.004 

044 ± 3 + 006 + 001 

U.U4+1 3= U.UUO _ .005 -0.005 

n no fi ,n 006 +a01 ° + 002 

U.UOU 33 U.UUU _o.011 -0.000 



Table 5: The measured gap fraction f (Arj) for the region x® BS < 0.75 unfolded 
with the average correction factors o/Pythia and Herwig. The statistical error, 
systematic errors, and calorimeter energy-scale uncertainty on the measurement 
are also listed. 
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W bin ( GeV) 


££UT GeV 


/ ± stat ± sys ± cal 


150.0,180.0 
180.0,210.0 

01 n n o/in n 
zlU.U, /4U.U 

240.0,260.0 


0.5 


0.019 ± 0.008 t° T 8 ±8Z 
0.013 ± 0.002 %™ $™ 

nniR _u n nno +0.006 +0.000 
0.016 ± 0.002 _ 0004 _ 0002 

0.021 ± 0.002 +g;g» 


150.0,180.0 
180.0,210.0 
210.0,240.0 
240.0,260.0 


1.0 


0.032 ±0.009 ±8;81 ±8:882 

027 ± 03 +0 - 004 +0 - 001 

U.UZf =C U.UUO _ .005 -0.002 

0.027 ±0.002 ±8:88! ±8:8S 
0.028 ±0.002 ±8:881 ±8:88f 


150.0,180.0 
180.0,210.0 
210.0,240.0 
240.0,260.0 


1.5 


nn7 7 + 014 +0 - 068 +aoo ° 

U.UM I U.U11 -0.058 -0.024 

n 044 ± 04 +0 - 005 +0 - 004 

U.U44: IC U.UU^i -0.005 -0.001 

04S ± 02 +0 - 006 +0 - 003 

U.U^O 3= U.UUZ _ .007 -0.004 

04S ± 02 +0 - 005 +0 - 002 

U.U^O 33 U.UUZ _ .005 -0.002 


150.0,180.0 
180.0,210.0 
210.0,240.0 
240.0,260.0 


2.0 


nii-i + nni5 +0 - 048 +0 - 000 

U.IIO 33 U.UIO -0.048 -0.018 

0.067 ± 0.004 +8:8o 3 r ±8Z 

0.069 ±0.003 +8Z 7 ^8:S§f 
0.064 ±0.003 ±8:881 ±:S:S8I 



Table 6: The measured gap fraction f (W) for the region x® BS < 0.75 unfolded 
with the average correction factors o/Pythia and Herwig. The statistical error, 
systematic errors, and calorimeter energy- scale uncertainty on the measurement 
are also listed. 
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Figure 1: The inclusive dijet cross section, differential in Etf Ap . The black circles 
represent the ZEUS data, with the inner error bars representing the statistical errors 
and the outer error bars representing the statistical and systematic uncertainties 
added in quadrature. The solid black line shows the prediction of Herwig and the 
black dashed line shows the prediction of Herwig plus BFKL Pomeron exchange. 
The dot-dashed line shows the prediction of Pythia and the dotted line shows the 
prediction of Pythia plus high-t photon exchange. The band shows the calorimeter 
energy-scale uncertainty. 
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Figure 2: The top plot is the inclusive dijet cross section differential in At], the 
middle plot is the gap cross section differential in Ar] requiring that E^ AF < 1 GeV , 
and the bottom plot is the gap fraction, f , as a function of Arj. Other details as in 
Fig.m 
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Figure 3: The gap fraction, f , as a function of A77 for different requirements 



on E^ AP . Other details as in Fig. Ql 
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Figure 4: The gap fraction, f , as a function of An for different requirements 
on E^ AP . The black circles represent the ZEUS data, with the inner error bars 
representing the statistical errors and the outer error bars representing the statistical 
and systematic uncertainties added in quadrature. The solid black line shows the 
prediction of Herwig plus BFKL Pomeron exchange. The open squares represent 
the HI data [5] scaled for comparison to the ZEUS phase space as described in the 
text. 
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Figure 5: The gap fraction, f, as a function of At] for different requirements on 
Etf AP . The black circles represent the ZEUS data, with the error bars representing 
the statistical and systematic uncertainties added in quadrature. The resummed 
calculation [6] is shown by the solid curve and the renormalization scale uncertainty 
is shown by the shaded band. The data are plotted at the 4 bin centers in Arj and the 
theory curve was produced by joining the bin centers for the ratios of the integrated 
cross sections for 8 bins in Arj. 
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Figure 6: The top plot is the inclusive dijet cross section differential in x 
the middle plot is the gap cross section differential in i^ BS requiring that E^ AP < 
1 GeV , and the bottom plot is the gap fraction, f , as a function of x® BS . Other 
details as in Fig. \J\ 
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Figure 7: The gap fraction, f , as a function of x° BS for different requirements 
on E^ AP . Other details as in Fig. [JJ 
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Figure 8: The top plot is the inclusive dijet cross section differential in W , the 
middle plot is the gap cross section differential as a function of W requiring that 
E^ AP < 1 GeV , and the bottom plot is the gap fraction, f , as a function of W . 
Other details as in Fig. d 



23 



ZEUS 



0.05- 







£ o.2|- e t ap < 5 GeV 

^ " • ZEUS (39 pb 1 ) 

HERWIG 

n 1t?L ■■■■ HERWIG + BFKLx 1.02 
U,igr . PYTHIA 

■■ PYTHIA + high-tyx 327 

■I 



In ■ ■ t^m n n n n A 



1 50 200 250 
W (GeV) 




1 50 200 250 

W (GeV) 



5 o.2- E? AP <1.0GeV 
**— 

0.15 
0.1 
0.05 




1 50 200 250 
W (GeV) 



0.2 
0.1 



I 


E? AP < 


2.0 GeV 


......... 



















1 50 200 250 
W (GeV) 



Figure 9: The gap fraction, f , as a function ofW for different requirements on 
E^ AP . Other details as in Fig. Ql 
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Figure 10: The top plot is the inclusive dijet cross section for x® BS < 0.75 
differential in Ar], the middle plot is the corresponding gap cross section differential 
in Ar/ requiring that E^ AF < 1 GeV , and the bottom plot is the gap fraction, f , as 
a function of At]. Other details as in Fig. Ui 
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Figure 11: The gap fraction, f , as a function of An for < 0.75 and different 
requirements on E^ AP . Other details as in Fig. [II 
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Figure 12: 27ie top plot is the inclusive dijet cross section for x® BS < 0.75 
differential in W , the middle plot is the corresponding gap cross section differential 
in W requiring that Etf AF < 1 GeV , and the bottom plot is the gap fraction, f , as 
a function ofW. Other details as in Fig. Ul 



27 



4.08 



0.04- 
0.02-4 



E? AP < 0.5 GeV 

n nc * ZEUS ( 39 P b 1 ) 
0.06I HERWlG 

HERWlG + BFKLx1.02 

... PYTHIA 

PYTHIA + high-tyx 327 



ZEUS 

^0.08 




150 



200 250 
W (GeV) 



-0.15- 



0.1- 



0.05 ■■ 



E? AP < 1 .5 GeV 



1 50 200 250 
W (GeV) 




1 50 200 250 

W (GeV) 



-0.15- 



0.1- 
0.05- 



E° AP < 2.0 GeV 



x° BS < 0.75 



1 50 200 250 
W (GeV) 



Figure 13: The gap fraction, f , as a function ofW, forx^ < 0.75 and different 
requirements on E^ AP . Other details as in Fig. [H 
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